Comment on "High Precision Measurement of the Thermal Exponent for the 
three-dimensional XY Universality Class" 
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A recent paper [Burovski et al, cond-mat/0507352] reports on a new, high-accuracy simulation 
of the classical <p i model (in the three-dimensional XY universality class). The authors claim that 
a careful scaling analysis of their data gives v = 0.6711(1) for the thermal critical exponent. If 
correct, this would neatly resolve the discrepancy between numerical simulations and experiments 
on 4 He. There is reason, however, to doubt the accuracy of the result. A re-analysis of the data 
yields a significantly higher value of v, one that is consistent with other Monte Carlo studies. 



Universality is an elegant and powerful concept, but as 
the basis for any method of data analysis it offers many 
dangers to the practitioner. For example, it is notori- 
ously difficult to extract critical exponents from data on 
finite systems. Such an analysis is exquisitely sensitive 
to the choice of scaling form, to the number of fitting pa- 
rameters, and to the range of system sizes included in the 
fit. Worse, it is quite difficult to quantify the uncertain- 
ties associated with these factors, so there is a tendency 
to overstate the accuracy of measured exponents. Refer- 
ence Q may suffer from this very problem. 

After repeated re-analysis of their experiments on the 
superfluid transition of 4 He in microgravity 0, 0, Ej], 
Lipa et al. have concluded that the best experimental 
value of the thermal critical exponent is v c ^ v — 0.6709(1). 
This value is not in good agreement with numerical 
simulations of other models also believed to be in the 
three-dimensional XY universality class. Two recent 
Monte Carlo studies dEI found u mc = 0.6723(3) (8) and 
0.67155(27). 

In an effort to settle the controversy, Burovski et al. 
studied the classical </> 4 model 
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by expanding the partition function in series and sam- 
pling the various terms stochastically. They considered 
two critical points — 

point A : t = -0.0795548(1), U = 0.4101562(14) 
point B : t = -0.07142883(7), U = 0.3605750(8) 

— the second chosen so as to minimize the subleading 
corrections to finite-size scaling (although one wonders if 
this helps or simply makes it harder to pick out the sub- 
leading corrections when it comes time to fit the data). 

Burovski et al. computed the superfluid stiffness p s and 
two of its t derivatives for a range of linear systems sizes 
L = 4, . . . , 96. Their results for R' = d{p s L)/dt, evalu- 
ated at the two critical points, are reproduced in TableQ] 
Each of the entries corresponds to 5 x 10 8 Monte Carlo 
sweeps. 



TABLE I: Derivative of the superfluid stiffness scaling func- 
tion at critical points A and B computed via Monte Carlo. 





R' = d{p s L)/dt\t =tc 


L 


data set A 


data set B 


4 


2.0329(9) 


1.9907(3) 




5 


2.8414(5) 


2.7843(7) 




6 


3.7316(4) 


3.6586(1)" 




7 


4.6955(5) 


4.6064(5) 




8 


5.7289(4) 


5.6221(9) 




9 


6.8265(7) 






10 


7.9848(9) 


7.840(1) 




11 


9.2031(13) 






12 


10.474(2) 


10.286(1) 




16 


16.074(4) 


15.789(2) 




20 


22.403(4) 


22.020(3) 




24 


29.396(7) 


28.897(3) 




32 


45.095(13) 


44.36(1) 




48 


82.48(3) 


81.15(2) 




64 




124.57(7)° 




96 


231.56(17) 







"points of disagreement between Fig. of this paper and Fig. 3 
of Ref.Q 



They simultaneously fit the two data sets, assuming a 
scaling form 



R' = (const.)£ 1/l/ [1 + (const. )L^] . 



(2) 



The constants were allowed to take on different values 
at critical points A and B, whereas the exponents v and 
uj were assumed to be universal (i/ = vx = vb, etc.). 
On the basis of a known irrelevant exponent, the fit was 
constrained such that \lu — 0.795 < 0.03; data for system 
sizes below £ C utoff = 12 were discarded. Final fit values 
of w — 0.796(3) and v = 0.6711(1) were reported. 

Is this value of v convincing? Such a high degree of 
confidence in the fourth digit is probably not warranted. 
It is clear even to the eye — in Figs. 2 and 3 of Ref.Q] — that 
the data allows considerable leeway in the slope and offset 
of the fit. Moreover, two aspects of the data analysis are 
troubling: 

1. It is not particularly useful to restrict the value 
of the scaling exponent. While it is true that the 
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FIG. 1: (Top panelVThe data from Tableware plotted exactly 
as m Fig. 3 of Ref. Q The two black squares mark the data 
points in Fig. 3 (at L = 6 and L — 64) that do not coincide 
with those shown here. The solid lines denote the best fits 
for critical points A (red) and B (pink) with L cuto ff = 12 and 
(j = 0.796 fixed. Best fits for L cuto a = 4,5,6,7,8,9,10,11 
are drawn dotted and those for L cuto g = 16, 20, 24 are drawn 
dot-dashed. (Bottom panel) The optimal v — computed with 
both u) held fixed and ui left free to vary — is plotted as a 
function of the lower size cutoff. The corresponding values of 
ui are shown in the figure inset. The violet band indicates one 
standard deviation above and below u exp from Ref. 0- The 
Monte Carlo results from Refs. |f| and Ref. are marked by 
the cyan and hatched regions. 




lytic in L 1 [which are not explicitly included in 
Eq. (J2J]. In practice, cu, as it appears in Eqs. (4) 
and (6) of Ref. Q, is an effective exponent that 
approximates the sublcading behaviour over some 
range of L. pfl 

2. More important, the authors have failed to quantify 
the effect of L C utoff on their fit. If, as I believe is 
true in this case, the optimal v depends sensitively 
on the choice of the lower size cutoff, then some 
convincing criterion must be advanced to justify 
choosing one value of £ C utoff over another. 

In the top panel of Fig. I have replotted the super- 
fluid stiffness data of Burovski et al. alongside my own 
best fits. These fits were generated using exactly the 
procedure outlined in Ref. yj for a range of L cuto ff , with 
and without constraints on w. The large variability in v 1 
as seen in the bottom panel of Fig. ^ suggests that an 
uncertainty of ±0.0001 is overly optimistic. Moreover, 
contrary to the authors primary claim, it appears that 
(for 6 < £ cu toff S 20) v is largely consistent with the 
result of Campostrini et al. and not significantly closer 
to the experimental value. There is little compelling ev- 
idence for a value as low as 0.6711. 

It is important to note that there is a curious disagree- 
ment (in two data points from the B data set ) between 
Fig. □ of this paper and Fig. 3 of Ref. 1. In Fig. 3, the 
L = 6 data point has the same position, but a larger 
crrorbar; the L — 64 data point is placed considerably 
higher. The origin of this discrepancy is somewhat mys- 
terious. It may simply be a consequence of minor ty- 
pographical errors in Table ID In that case, Burovski et 
al. are fitting to slightly different data set than I am. 
Nonetheless, it is telling that minor changes on the order 
of one or two standard deviations to two of twenty-eight 
data points could have such a disruptive effect. It sug- 
gests that v = 0.6711(1) is an unrealistic estimate of the 
thermal critical exponent. 



lowest-order subleading correction arises from an ir- 
relevant scaling field with exponent uj- m ~ 0.8, that 
correction coexists with other contributions ana- 



[1] 
[2] 
[3] 



[5] 

[6] 

[7] 



E. Burovski, J. Machta, N. Prokofev, and B. Svistunov, 
cond-mat/0507352. 

J. A. Lipa, D. R. Swanson, J. A. Nissen, T. C. P. Chui, 
and U. E. Israelsson, Phys. Rev. Lett. 76, 944 (1996). 
J. A. Lipa, D. R. Swanson, J. A. Nissen, Z. K. Ceng, P. R. 
Williamson, D. A. Strieker, T. C. P. Chui, U. E. Israelsson, 
and M. Larson, Phys. Rev. Lett. 84, 4894 (2000). 
J. A. Lipa, J. A. Nissen, D. A. Strieker, D. R. Swanson, 
and T. C. P. Chui, Phys. Rev. B 68, 174518 (2003). 
M. Hasenbusch, T. TorSk, J. Phys. A 32, 6361 (1999) 
M. Campostrini, M. Basenbusch, A. Pelissetto, P. Rossi, 
and E. Vicari, Phys. Rev. B 63, 214503 (2001). 
K. S. D. Beach, L. Wang, and A. W. Sandvik, 
cond-mat/0505194. 



